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Abstract 

We propose the general scaling model for the diffusion- annihilation reaction 
Af + A- — ► with long-range power-law interactions. The presented scal- 
ing arguments lead to the finding of three different regimes, depending on 
the space dimensionality d and the long-range force power exponent n. The 
obtained kinetic phase diagram agrees well with existing simulation data and 

approximate theoretical results. 
PACS: 82.40, 61.30J 
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The problem of annihilation-diffusion kinetics has generated significant interest in recent 
years, since it was shown [T|] - || that in a bimolecular reaction A + B — > 0, long- wavelength 
density fluctuations cause significant slowing down of the particle density decay. While the 
ordinary kinetic-rate approach suggests that in a bimolecular reaction, particle density p 
always decays as t' 1 , the account of fluctuations leads to a power-law decay: 

p<xt- v , (1) 

where u < 1; when no long-range interaction exists between particles, v = d/4 for 
d < 4 and v — 1 when d > 4, where d is space dimensionality [Q - @. Thus, 4 is critical 
dimension for the bimolecular reaction, and below it, long- wavelength fluctuations determine 
annihilation kinetics entirely. 

Even though the role of fluctuations in systems without long-range interactions is well 
understood and confirmed by numerical simulations 0-0, the situation is far less clear 
when long-range forces are present. Meanwhile, the latter problem is especially important 
in different physical applications, e.g., in condensed matter physics, when one replaces par- 
ticles A and B with dislocations or vortices and analyzes their annihilation upon quench 
from a high-temperature to a low-temperature phase. These considerations prompted sev- 
eral numerical studies into the problem of annihilation kinetics in Coulombic systems in two 
dimensions f|, @, [fLOR . While it was shown that power-law decay ([!]) is satisfied on a late 
stage of the annihilation process, and that exponent v is rather close to 1, its value has not 
been determined fully or calculated rigorously. In our recent work ||, numerical simulations 
in two dimensions yielded v = 0.9 ± 0.05. The scaling theory was proposed to explain this 
unusual result ||, JTTJ] , based on an assumption that the system first relaxes charge density 
fluctuations towards statistical equlibrium and only then proceeds with annihilation. How- 



ever, other theories [|12[] - |13[, predicted v = 1 for Coulombic systems, if initial conditions 
remain random when annihilation starts. 

Quite recently, we proposed a self-consistent description of the annihilation-diffusion 
kinetics p| , which allows to take into account different initial conditions. This theory 



neglects fluctuations of the total particle density, but, unlike a mean- field theory, takes into 
account the fluctuations of the charge density Such an approximation enables one to solve 
equations of evolution and determine asymptotics for different scaling regimes. In the case 
of Gaussian random initial charge density, kinetic phase diagram was calculated and final 
(large-t) asymptotics were as follows: 

if n > 1 + d/2 and d < 4 (fluctuation-dominated region) 

p oc (Dt)-", v = d/4; (2) 

if n < 1 + d/2 and n > d — 1 (intermediate region) 

P (x(Qt)- u , v= f ; (3) 

H vv ; ' 2-d + 2n K ' 

and if d > 4 (mean-field region) 

p oc (£t)-", v = 1, (4) 

with D, Q and JC being diffusion, electrostatic and annihilation constants. Equations (^) 
- (^) do not exhaust all scaling solutions of self-consistent equations of evolution, these 
are only the large-t asymptotics. Equations (0) and (|J) are very well-known, while the 
intermediate solution @ was proposed in the scaling theory of Ispolatov and Krapivsky [13| 
(but the region of its applicability was determined in a different manner). 

In this Letter, we propose a simple scaling theory to derive relations (fj) - (f|). Such a 



theory would not only support the validity of the self-consistent approach of Ref. ||14|| , but 
also would give an additional insight on the role of initial conditions in determining late 
stage asymptotics for annihilation-diffusion problems. 

Let us start by writing a Langevin equation for the i-th particle, assuming that inertia 
forces are negligible in comparison with viscous drag force: 

dri ^ Qqiqj 

^ = ^K^F T(ri " rj) + fri ' (5) 
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where (3 is a "friction coefficient", and f r i is a random Brownian force acting on the i-th 
particle. We will use the Langevin equation to characterize a slow relaxation of an initial 
particle density fluctuation and thus estimate the annihilation rate. 

Let us select a region with characteristic size L. Because of the assumption about ran- 
domness of the initial particle distribution, 5N oc y/N, where N is the number of particles 
inside the volume V = L d . The mean-square normalized density fluctuation in this volume 
is: 

5 -P oc L- d '\ (6) 
Po 

As it is traditionally done in scaling theories of annihilation, we assume that excess 
particles start annihilating only when they travel the distance on the order of L. In this 
case, assuming that diffusion is irrelevant, we rewrite the Langevin equation, replacing all 
occurences of r with L: 

dL _ 5_Q(L)_ 

1 dt L n ' K) 

where total charge fluctuation in the region SQ(L) oc L d//2 , again due to the Gaussian 
initial conditions. From this follows: 

L n+l-d/2 K ^ (g) 



or 



Loctl +1 - d/2 , (9) 

where ti represents time necessary for an average particle to travel the distance L (as- 
suming it does not annihilate during this time). At the time t^, majority of particles that 
originally were in the volume V, have annihilated, and only those excess particles still re- 
mained, so the density is: 

PM « 5 J1 oc L -d/2 oc t ~^T^m^ (1Q) 
Po A 



or, omitting the index " L" , 

P OC t~~ 2n + 2-d . (11) 
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In order for scaling laws ([8]) and ( |iT|) to be valid, it is necessary that i) drift occurs slower 
than annihilation, and that ii) drift occurs faster than thermal (Brownian) diffusion. These 
conditions imply that: 

n + l-d/2<2, (12) 



11 < 1, (13) 



2n + 2 - d 

which is equivalent to: d — 1 < n < 1 + d/2. 

If condition (|T2"D is not satisfied, i.e., n > 1 + d/2, deterministic drift is irrelevant, and 
everything is determined by the competition of diffusion and annihilation. In this case, as 
expected: 

L oc t 1 / 2 , (14) 
p oc r d/ \ d<4; (15) 

pocr 1 , d>4. (16) 

Thus, using simple scaling approximations, we confirmed validity of the self-consistent 
approximation for the annihilation-diffusion system with long-range interactions. Scaling 
exponents obtained in the large-t limit agree well with known theoretical models and nu- 
merical simulations. We also showed that the initial density distribution plays an important 
role in determining large-t scaling behavior at least for strongly interacting systems such as 
Coulombic, and made estimates of this behavior when the initial distribution is Gaussian. 

In the future, more theoretical and numerical efforts will be required to explore in more 
details the rich kinetic phase diagram for systems with long-range interactions. While the 



use of scaling arguments like the one proposed here cannot be used for exact analysis of 
annihilation behavior, it can and will complement other numerical and analytical methods 
in elucidating main features and dependencies of annihilation kinetics. 

Acknowledgement. This work was done during the March Meeting of the American 
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